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**F ILE**1D**MTHGCOSH 
MM MM TTTTTTTTTT HH HH GGGGGGGG cccccccc 000000 SSSSSSSS HH HH 
MM MM TTTTTTTTTT HH HH GGGGGGGG cccccccc 000000 SSSSSSSS_ HH HH 
MMMM = MMMM TT HH HH GG cc 00 00 SS HH HH 
MMMM MMMM TT HH HH GG cc 00 00 SS HH HH 
4M 6M UMM TT HH HH GG cc 00 00 SS HH HH 
MM MM UMM TT HH HH GG cc 00 00 SS HH HH 
MM MM TT HHHHHHHHHH = GG cc 00 00 SSSSSS HHHHHHHHHH 
MM MM TT HHHHHHHHHH =GG cc 00 00 SSSSSS HHHHHHHHHH 
MM MM TT HH HH GG GGGGGG CC 00 00 SS HH HH 
MM MM TT HH HH GG GGGGGG CC 00 00 SS HH HH 
MM MM TT HH HH GG GG CC 00 00 SS HH HH coos 
MM MM TT HH HH cc 00 00 SS HH HH cece 
MM MM TT HH HH GGGGGG cccccccc 000000 SSSSSSSS HH HH eeee 
MM MM TT HH HH GGGGGG cccccccc 000000 SSSSSSSS HH HH cove 
LL IIIT] SSSSSSSS 
LL HII! SSSSSSSS 
LL I] SS 
LL I] SS 
LL I] $$ 
LL I] SS 
LL I] SSSSSS 
LL II SSSSSS 
LL I] SS 
LL I] SS 
LL I] SS 
LL I] Ss 
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; if paoaring Hyperbolic Cosine routine 
: File: MTHGCOSH.MAR EDIT: RNH1005 


-TITLE MTHSGCOSH 
- IDENT /1-005/ 


COPYRIGHT (c) 1978, 1980, 1982, 1984 B 
DIGITAL EQUIPMENT _GORPORATION, MAYNARD. MASSACHUSETTS. 
ALL RIGHTS RESERVE 


OPYRIGHT NOTICE. THIS SOFTWARE OR ANY OTH ER 
E PROVIDED OR OTHERWISE MADE AVAILABLE TO ANY 
E TO AND OWNERSHIP OF THE SOFTWARE IS HEREBY 


H 

THE TERMS OF SUCH LICENSE AND WITH THE 
B 

OTHER PERSON. NO TI 
TRANSFERRED. 

THE INFORMATION IN THIS SOFTWARE IS SUBJECT TO CHANGE WITHOUT NOTICE 
ae paket ite NOT BE CONSTRUED AS A COMMITMENT BY DIGITAL EQUIPMENT 


DIGITAL ASSUMES NO RESPONSIBILITY FOR THE USE OR RELIABILITY OF ITS 
SOFTWARE ON EQUIPMENT WHICH IS NOT SUPPLIED BY DIGITAL. 


REAR AAAE AA AAAAAAAARARAARAAARAAAEAAAAEAAAAEAAERAEAAAAEAEAAAEAAEAEEAAAEEEE 


FACILITY: MATH LIBRARY 
ABSTRACT: 
MTH$GCOSH is a function which returns the G eae ben hyperbolic cosine 


of its G floating point argument. The call is standar 
call-by-reterence. 


VERSION: 1 


HISTORY: 
AUTHOR: 


Steven B. Lionel, 26-Jan-79: Version 1 


MODIFIED BY: 
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-SBTTL HISTORY ; Detailed Current Edit History 
dit History for Version 1 of MTH$GCOSH 


1 = Adapted iros gag” oa 
Use MTHSGEXP_R6. i750 £197 
B aven'y sécond 7 . ¢ Ie say » neyt —- between 25 and 

tln2 by computing COSHex) = (Z + = EXPCix!) 
Et taioste ed econd call to EXP for Rd. values between 27#Ln2 
subrira att final floating point Ohytetons by 2 to interger 
subrtracts of 1 from the ttgegs RP = 

Changed “nat mask to exclu yee R7 = pb longer needed. 
oe = mum range to 1 
chanege ed Lo ic for computing exbeiat~in2) to reduce error. 
Changed oe to G* in call to MTH$$SIGNAL RNH 09-Sept-1981 
Eliminated symbolic short literals. RNH 15-Oct-8 
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, -SBTTL DECLARATIONS 3; Declarative Part of Module 
f ; INCLUDE FILES: 
fh > EXTERNAL SYMBOLS: 
6 *EXTRN MTHSSSIGNA 
7 eEXTRN MTHSK ELOOVEMAT 
4 SEXTRN MTHSGEXP_ 
0 : EQUATED SYMBOLS: 
00000004 § value =4 3; value.rg.r 
00 4; 
00 5 ; MACROS: none 
4 $ 3 
0000 BB 
0000 9; 
4 $ 3 PSECT DECLARATIONS: 
veag 3 -PSECT _MTHSCODE PIC,SHR,LONG, EXE ,NOWRT 
4444 as 3 program section for math routines 
0900 95 OWN STORAGE: none 
00 96 ; : 
44 97 ; 
00 98 : CONSTANTS: 
44 99 ; 
000 100 
000 101 G6_1023_L0G_2: 
7B606E3A 288740A6 000 1 § -QUAD “%*X7B606E3A28B740A6 3 1023*ln2 
3 103 6_1024_L0G_2: 
39EEFEFA 2E4240A6 , I 6.LOG.2. ar “X39EEFEFAZE4240A6 > 1024*Ln2 
3COOFEFA 2E424006 a I 6_LOG.2 toe “X3COOFEFA2E424006 3 Chigh 43 bits of Ln2)+2**#-43 
4C67361C 86548050 001 1 8 Re “%4(67361C86548D50 : ln2 = G_LOG_2_HI 
109 G_2_ POWER 7: 
00000000 000041C0 119 UAD %*x00000000000041C0 3 2e*27 
11 GCOSHTAB: 
850A 3682 11 .WORD “O08 02010881 
7A 67DD f 114 ~ WORD 378 "011557 3: DECIMAL: 0.2759648863787355D-06 
O19E SFIA 115 WORD 72000083 
194A 3565 4 138 WORD “0 2 "001451 3 DECIMAL: 0.2480155975461668D-04 
16C 3F 8 11 - WORD 701393 
a 7 f 18 -WORD *0 13400,°0 ; 1 3; DECIMAL: 0.1388888889781712D-02 
555 aS 4 11 ~ WORD H 705,°0 ¢3 5 
4DF9 55 44 120 -WORD * §25,°0046771 3 DECIMAL: 0.4166666666665359D-01 
4 48 121 eWORD *004 9 
1 4 1 ¢ ~WORD Q “0 6 3 DECIMAL: 0.500000000000000000 
4 1 WORD 694 0 
124 -WORD 0, 3 DECIMAL: 0.1000000000000000D+01 
00006 8 125 GCOSHLEN = .- GCOSHTAB/8 
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-SBTTL MTHSGCOSH = G Double Precision Floating GCOSH 


~;| 


+4 

; FUNCTIONAL DESCRIPTION: 

: GCOSH = double precision G floating point function 
; GCOSH(X) is computed as: 


If ixi < 2e-27 x) = 1, 


: GCOSH( 

; IF 2ee-27 =< Ix! <0. 23. GcOsiicx) = Chebyshev series 

: If 25 < Ix! < 37 sing t 2 = GEXP(IX!) and set GCQSH(X) = (241/2)/2 

: is n2-s< Ix! 1023ein2,. then GCOSH(X) = GEXPUIKID Jo 

; tin2 =< IX! < 1024*Ln2, then GCOSH(X) = GEXP(!X!-Ln2). 

3 if 19 4*ln2 =< IX!i, then overflow. 

: CALLING SEQUENCE: 

: GCOSH.wg.v = MTHSGCOSH(x.rg.r) 

: INPUT PARAMETERS: 
00000004 LONG = 4 3; define longword multiplier 
00000004 x = 1 * LONG S Contents of x 18 the ergunent 


IMPLICIT INPUTS: none 
OUTPUT PARAMETERS: 
VALUE: G floating hyperbolic cosine of the argument 
IMPLICIT OUTPUTS: none 
COMPLETION CODES: none 
SIDE EFFECTS: 
Signal: MTHS_FLOOVEMAT if 1024*Ln2 c xt gy reserved operand in RO/R 
(copied to the signal mechanism vector CHFSL Mgt i byt LIBSSIGNAL). 
Associated message is: ‘FLOATING OVERFLOW IN “RATA Result is reserved 
operand -0.0 unless a user supplied (or any) error ’eanalae “changes CHFSL_MCH_RO/R1 


wri re procedure disables floating point underflow, enables integer 
ver 


407C ENTRY MTHSGCOSH, “M<IV, R2, R3, R4, RS, RO > 
3 standard call-by-reference entry 
:; disable DV (and FU), le IV 


MTHSFLAG_JACKET 3 flag that this is a jacket procedure in 
MOVAB G*MTHSSJACKET_HND, (FP) 


oe et a a 4 — — 4) —) —4 9 —) —4 —9 —4 —_» - 2s — 2 2 2) —) 2) —) ) — 9 2 os 1 a 3 1 ss a a 


BAIN PAAAAAA AAD MMIII PPP PLL EE 


SOON UE WIN SO OD NAWULE WN 9 OD NAN EWN OOD NAU EW OOO NAUE W000 


6D 00000000°GF 9€ 


set handler address to jacket 
handler 
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50 Ss C SOFD 
HY g0 0 BF AA 
FFO 8F 9 1 
1 8 

3€60 8F 50 ©B1 
0S = 18 


50 08 SOFD 
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oul 


FF71 CF 50 5S1FD 
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GEQ_T0_2M27: 
MULG2 

POLYG 

RET 
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Xi 
GEQ_T0_6.25: 
CMPG 
BGTR 


ONE_TERM_ONLY: 
UBW 
RET 


3 2ee-27 =< 1X! 


1023*ln2 < ix! 


C 11 


c Cosine gehen 16-SEP 
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RO otanse 


GEG_TO_2™ 


#1, RO 


< 0.25 


RO,RO 
RO, #GCOSHLEN~1 


G_1023 


RO 062 
GTR_ TRAN. 1023, C0G_2 


; 0.25 =< IX! =< 1023#Ln2 


MTHSGEXP_R6 


0, G2 POWER _27 
aia 


#*x0010, RO 
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case of an error in routine 
If an error, convert signal to user PC 


+ and resigna 

3 RO/R1 = (Xi = @value (AP) 
3 RO/R1 = IX! 

3 compare {Xi with 8-32 

: branch if {xi >= 0.25 

3 compare {Xi with iat 
3 branch if IX! >= 2e%-27 
: RO/R1 = 1.0 


return with result = 1.0 


;Get ARG**2 for POLYG. 


3; RO/R1 = SUM(Ci*X**7) 
: return with result in RO 


: compare {Xi with 1023*Ln2 
: branch if {Xi > 1025*Ln2 


RO/R1 = GEXPC; 

Compare GEXP(} with 2**27, if 
Larger. onty one term is needed. 
R2/R3 = GEXP(- a 


RO/R1 = GEXP(X EXP (=X) 


RO/R1 = (GEXP(X) + GEXP(-xX))/2 
return with result in RO/R1 


Xi) 
Xi) 
ne 

Xi 
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GTR_THAN_ 1023_L0G_2: 
CMPG ~ RO 


-1024_L06_2 3; compare {Xi with 1024*\n2 
BGEG ERROR 


Sel 
branch to ERROR if 1 $24 eln2 =<ix! 


1023*ln2 =< ix! < 1024*Ln2 


SUBG2 G_LOG.2_HI ; RO/RI=IX!-(high 43 bits of ln2+2e*-43) 
JSB MTHS GEXP re * ; RO/RI=GEXP(iXT=G_LOG 2 HI) 

MULG3 tse 2 Co: RO, R2 3 R2/R3=G_LOG_2 LOSGEXPCT TX!=G_LOG_2_LO) 
ADDG2 Rot RO ~ 3 RO/R1=GEXP(TXT-Ln2) 

RET : return with result in RO/R1 


$ 1023*ln2 + LOG(2) =< 1X1, error 


ERROR: MOVZBL #MTHSK_FLOOVEMAT, =-(SP) 
ASHQ #15, #1, R 


condition value 
RO = tags = reserved operand -0.0 
goes +s spe\ wer vector 
CheSL OM MCH + ag ) so error handler 
cen modify the result. 
_— error and use real user's PC 
in ependent of CALL vs J 
return - RO/R1 restored from CHFSL_MCH_RO/ 


CALLS #1, G*MTHSSSIGNAL 
RET 
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ERROR R 01 
GCOSHLEN = 
GCOSHT R 1 
GEQ_T0_0. ; R 1 
CTRTINAN, 1023_L0G_2 BR Oy 
ee 1 
G_2_POGER_2 R 1 
G_LOG_2_HT R 1 
G LOG_2_L0 18 R 1 
LONG = 04 
MTHSSJACKET_HND eeeeeeee =X 01 
MTHSS$SIGNAL rrereeee x 0 
MTHSGCOSH 00000058 RG 1 
MTHSGEXP_R6 keeerene K 0 
MTHSK_FLOOVEMAT teeeeeee =X 0 
ONE _TERM_ONLY 000000A7 R 1 
VALUE = 00000004 

posemeosemaesneeoaes + 

; Psect synopsis ! 
PSECT name Allocation PSECT No. Attributes 
» oe. 00000000 «¢ _ 0.) 00 ¢ 0.) NOPIC USR CON ABS LCL NOSHR NOEXE NORD NOWRT NOVEC BYTE 
_MTHSCODE 000000DB ¢ 219.) 01¢ 1.) PIC USR CON REL LCL SHR EXE RD NOWRT NOVEC LONG 

Sewcceousesscesuessesones} 
: Performance indicators ! 

Phase Page faults CPU Time Elapsed Time 
Initialization 29 0:00:00.0 00:00:00.73 
Command processing 114 8: : 8 :00:04.31 
Pass 1 89 :00:00. : 3: 3.70 
Symbol table sort 8 0:00: 3} :00:00.01 
Pass 2 6 + :00.71 - 8: 4-38 
Symbol table output 4 :00: 83 :00:00.0 
Psect synopsis output :00:00. 0: 8: 04 
Cross-reference output :00: 9 0:00: 90 
Assembler run totals 30 0 ‘ 0:00:10.9 


{he working set Limit was 900 pages. 
3992 bytes (8 pages) of virtual memory were used to buffer the intermediate code. 

There were 10 pages of symbol table space allocated to hold 20 non-local and 0 local symbols. 
source Lines were read in Pass 1, producing 11 object records in Pass 2. 

1 page of virtual memory was used to define 1 macro. 
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! Macro Library statistics ! 


$e erm ese ne new eee nem e eee ecw cen} 


Macro Library name Macros defined 
_$255$DUAZ8: [SYSLIBISTARLET.MLB; 2 Carmine 
0 GETS were required to define 0 macros. 

There were no errors, warnings or information messages. 


MACRO/ENABLE=SUPPRESSION/D1SABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:MTHGCOSH/OBJ=OBJ$:MTHGCOSH MSRC$:MTHJACKET/UPDATE=(ENHS$ : MTHJACKET) #+MSRC 


NT CORPORATION 
D PROPRIETARY 
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